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Abstract. Chain of calculations which have to be performed 
to predict any kind of signal in a deep underwater/ice neu- 
trino detector necessarily includes the lepton propagation thro- 
ugh thick layers of matter, as neutrino can be observed only 
by means of leptons (muons, first of all, due to their large 
ranges) that are generated in IN interactions. Thus, 

the muon propagation plays a key role when analyzing data 
and it is important to understand clearly how transportation 
part of simulation chain contributes to total inaccuracy of fi- 
nal results. Here we consider sources of uncertainties that ap- 
pear in Monte Carlo algorithms for simulation of muon trans- 
port. The trivial but effective test is proposed to measure the 
propagation algorithm accuracy. The test is applied to three 
MC muon transport codes (PROPMU, MUSIC, MUM) and 
results are reported. 



1 Introduction 

The main challenge for existing and planned underwater/ice 
neutrino telescopes (Barwick et ai, 1991), (Ahrens et ai, 
2001), (Aslanides et ai, 1999), (Carmona et ai, 2001), (So- 
kalski and Spiering, 1992), (Balkanov et ai, 2000), (Resva- 
nis etai, 1 994), (Grieder etai, 200 1 ) is detection of neutrino 
of extraterrestrial origin. But neutrino is neutral weakly in- 
teracting particle and one can not observe it directly but has 
to detect it by means of lepton that appears in v N IN 
interaction and passes a distance in medium before being de- 
tected. The case with / = /i^ is considered, first of all, be- 
cause muons possess the best ability to propagate large dis- 
tances in contrast to (due to their short radiation length) 
and (due to their short life time) in a wide range of ener- 
gies. On the other hand, atmospheric muons (that represent 
a penetrating component of atmospheric showers resulting 
from interaction of primary cosmic radiation with the Earth 
atmosphere) are the principal background for neutrino sig- 
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nal. To eliminate it one has to predict correctly the detec- 
tor response on atmospheric muon flux. Also one should 
keep in mind the fact that atmospheric muons are of scien- 
tific interest themselves: for instance, from the point of view 
of charm production (Sinegovskaya and Sinegovsky, 2001). 
The last (but not the least) important point is detector cali- 
bration which has to be done, again by means of atmospheric 
muons that are the only more or less known intensive cali- 
bration source for deep underwater/ice neutrino telescopes. 

Thus, when analyzing experimental data obtained with un- 
derwater/ice detectors one has to apply some model for muon 
transport through thick layers of matter, using an algorithm 
that gives the muon energy Ex at the end of distance D pro- 
viding that its initial energy is equal to Ea . Commonly, Monte 
Carlo (MC) technique is used for this purpose (Aarnio et ai, 
1983), (Lipari and Stanev, 1991), (Antonioh et ai, 1997), 
(Kudryavtsev et al, 1999), (Sokalski et ai, 2001), (Rhode 
and Chirkin, 2001) being the most adequate to essentially 
stochastic nature of muon energy losses. Such a MC algo- 
rithm necessarily incorporates some uncertainties into the fi- 
nal result. These uncertainties can be divided into two parts; 

(A) insurmountable uncertainties which relates to 

- finite accuracy of formulae for muon cross-sections; 

- data on medium density and composition. 

(B) "inner" errors that are produced by simulation algorithm 
itself due to 

- finite accuracy of numerical procedures; 

- simplifications that are done to get reasonable com- 
putation time. 

In case {B) > (A) an algorithm is obviously too inaccurate 
and corrupts result more than uncertainties with the muon 
cross-sections and medium composition. (B) <C (A) is bet- 
ter but leads to unreasonable wasting of human and computer 
efforts, the tool for the muon transport is too "thin" compar- 
ing to insurmountable uncertainties of group (A). Thus, the 
case (B) < (A) seems to be an optimum equilibrium. 
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Fig. 1. Total cross-section ^7^* as a function ofW^ as measured and 
parametrized by ZEUS Collaboration (markers and dashed line, cor- 
respondingly), and parametrized by Bezrukov-Bugaev (solid line). 
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Fig. 2. Muon energy losses due to photonuclear interaction in pure 
water as obtained using (a) BB parameterization for cr^ot (solid line) 
and (b) ZEUS parameterization (dashed line). 



This report is aimed to get some quantitative parameters of 
muon transport algorithm which condition {B) < (A) leads 
to. 

2 How precisely do we know the muon enei^y losses? 

The accuracy of existing formulae for muon cross-sections 
was analyzed in number of works (KokouUn and Petrukhin, 
1991), (Rhode and Carloganu, 1998), (KokouUn, 1999). The 
conclusion is that it is not higher than wl% for i?^ < 1 TeV 
and becomes the worse the higher muon energy is. We do not 
give here the complete analysis but consider only two exam- 
ples that concern the muon photonuclear interaction which 
contributes the largest uncertainty to the total muon energy 
losses comparing with other kinds of muon interactions that 
represent purely electromagnetic (and hence, better known) 
processes: ionization, bremsstrahlung, and direct e+e~-pair 
production. 

For last twenty years one has been using the formula for 
photonuclear muon cross-section that was developed in the 
frame of generalized vector dominance model (GVDM) (Be- 
zrukov and Bugaev, 1981a). In such an approach the dif- 
ferential cross-section (and, consequently, the muon energy 
losses due to photonuclear interactions) is proportional to the 
total cross-section for absorption of a real photon by a nu- 
cleon, cr^jJJ. There are several parameterizations for a^^^ that 
were obtained by trying to get the best fit to experimental data 
(Bezrukov and Bugaev, 1981a), (Abramowicz et al., 1991), 
(Donnachie and Landshoff, 1992), (Breitweg et al., 1999), 
(Butkevich and Mikheyev, 2001). Two of them are shown in 
fig. 1 along with the last results from the ZEUS experiment 
(Breitweg et al., 1999). cr/oj is presented as a function of 



W^, where W is the center-of-mass energy of the -yp sys- 
tem. Parameterizations (Bezrukov and Bugaev, 1981a) (BB) 
and (Breitweg et al., 1999) (ZEUS) both fit more or less the 
experimental data but even within "experimentally known" 
range W'^ < 10^ GeV^ the difference between two param- 
eterizations reaches 4%. At W'^ = 10^ GeV^ the BB curve 
is higher than the ZEUS one by 15%. The same data are 
presented in fig. 2 in terms of muon energy losses due to pho- 
tonuclear interaction in pure water. Difference is kept within 
several percents below » 100 TeV and then grows up to 
w 10% at Ef^ = 1 EeV. Parameterizations for cr^j^ published 
in (Abramowicz et al., 1991), (Donnachie and Landshoff, 
1992) are lower comparing to ZEUS parameterization since 
they are based on old results which differ remarkably from 
the latest ZEUS data. On the other hand, the recent param- 
eterization (Butkevich and Mikheyev, 2001) is higher com- 
paring to BB. 

It became clear in last few years that at very high lepton 
energies the essential part of photonuclear cross-section is 
due to nonOVDM contribution. It appears, that this addi- 
tional part is well described by QCD perturbation theory. 
Recent new calculations (Dutta et al., 2000), (Butkevich 
and Mikheyev, 2001), (Bugaev and Shlepin, 2001) show an 
essential increase of the muon energy losses due to photonu- 
clear interactions. In (Dutta et al., 2000) the ALLM for- 
mulae (Abramowicz et al, 1991) (that is based on Regge 
approach and on HI and ZEUS data) is used for parame- 
terization of nucleon structure function F2. In (Butkevich 
and Mikheyev, 2001) the CKMT Regge model (Capella et 
al, 1994) is apphed for a description of F2 at low and inter- 
mediate Q"^ (squared four-momentum transfer) and the fit of 
parton distribution functions given by MRS group (Martin 
et al, 2000) at high Q^. Authors of (Bugaev and Shlepin, 
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Fig. 3. Muon energy losses in standard rock due to photonuclear 
interaction with- (Dutta et al, 2000), (Butkevich and Mikheyev, 
2001), (Bugaev and Shlepin, 2001) and without (Bezrukov and Bu- 
gaev, 1981a) accounting for QCD part (dashed lines and open cir- 
cles). Curves for bremsstrahlung (Bezrukov and Bugaev, 1981b), 
(Andreev et ah, 1994), (Kelner et al., 1995), direct e'^e^-pmr pro- 
duction (Kokoulin and Petrukhin, 1969), (KokouUn and Petruldiin, 
1971), (Kelner, 1997), (Kelner etal, 1999) (dotted lines) are shown, 
as well. Solid lines stand for total energy loss (brem + pair -I- nuc, 
ionization not included) with- (Dutta et at., 2000), (Bugaev and 
Shlepin, 2001) and without (Bezrukov and Bugaev, 1981a) ac- 
counting for QCD part in photonuclear interaction. 



2001) use for calculation of QCD perturbative part the color 
dipole model in version of (Forshaw et al., 1999), with pa- 
rameters founded from latest DESY data. 

The results of these three works are presented in fig. 3 and 
tables 1,2 in terms of the muon energy losses in standard 
rock. (Butkevich and Mikheyev, 2001) and (Bugaev and 
Shlepin, 2001) are close to each other while the muon en- 
ergy losses obtained in (Dutta et al, 2000) are essentially 
lower. We believe that this is due to the fact that old data 
(Abramowicz et al., 1991) were used in (Dutta et al., 2000) 
to fit the parameters of model. Anyway, the total muon en- 
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+ 10% 


+ 16 % 
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+ 29 % 


lOPeV 


+ 34% 


+ 42 % 




100 PcV 


+ 53 % 


+ 71 % 




1 EcV 


+ 81 % 


+ 119 % 





Table 1. Differences in terms of the muon energy losses in standard 
rock due to photonuclear interaction (in percents) between three 
works that account for QCD part (Dutta et at., 2000), (Butkevich 
and Mikheyev, 2001), (Bugaev and Shlepin, 2001) and Bezrukov- 
Bugaev formula based on GVDM (Bezrukov and Bugaev, 1981a). 





Dutta et al. 


Bugaev & Shlepin 


100 TeV 


+ 1 % 


+ 2% 


IPeV 


+ 3% 


+ 4% 


10 PeV 


+ 6% 


+ 7% 


100 PeV 


+ 11 % 


+ 14 % 


lEeV 


+ 18% 


+ 27% 



Table 2. Differences in terms of total muon energy losses in 
standard rock (in percents) between two works that account for 
QCD part (Dutta et at., 2000), (Bugaev and Shlepin, 2001) and 
Bezrukov-Bugaev formula based on GVDM (Bezrukov and Bug- 
aev, 1981a). 



ergy losses become higher when accounting for perturbative 
QCD part in photonuclear interactions by 1-^2% at i?^ = 100 
TeV and by 20^30% at = 1 EeV Fig.4 shows changes 
in muon survival probabilities in sea water and standard rock 
which incorporation of nonGVDM QCD-corrections results 
in ^. It is seen that for extremely high energies (E"^ > 1 PeV) 
effect is very significant. 
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Fig. 4. Muon survival probabilities vs slant depth in standard rock 

(upper panel) and sea water (lower panel) when simulating the 
muon propagation with photonuclear interaction treated by GVDM 
approach (Bezrukov and Bugaev, 1981a) (dashed lines) and when 
using also correction for QCD part as obtained by (Bugaev and 
Shlepin, 2001) (solid lines). Six pair of curves on each panel cor- 
respond to six initial muon energies from 1 TeV up to 1 EeV (from 
left to the right). 

Energy losses are proportional, in particular, to medium 
density. The density of water (ice) is known or, in any case, 

'Accounting for nonGVDM corrections becomes even more im- 
portant for r-leptons as photonuclear interaction dominates in r's 
energy losses at Er > 100 TeV 

^Data for this figure, as well as for fig.5-10, were obtained with 
the MUM code for /i and r-lepton transport (Sokalski et at., 2001) 
that is available by request to authors. 
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can be measured precisely. But since the main technique to 
select the neutrino induced events is looking for up-going 
muons from the lower hemisphere, one should know also 
the density and composition of sea bed rock. When this is 
unknown one has to use a standard rock model (p = 2.65 
g/cm"^, A = 22, and Z = 11). At the same time the measured 
rock densities at Gran Sasso and Frejus underground labo- 
ratories are equal, for instance, to pcranSasso = 2.70 g/cm^ 
and pFrejus = 2.74 g/cm^, respectively (Ambrosio et al, 
1995), (Rhode, 1993) which is 2-=-3% higher. Investigation 
of the Baikal lake bed (NT-200 neutrino telescope) showed 
a complex structure with several layers of different compo- 
sition and density that varies with increase the depth from 
PBaikai = 1-70 g/cm^ to pBaikal = 2.90 g/cm'"^ (±25% vari- 
ations) (Panfilov, 2000). Thus, error that is incorporated to 
muon transport calculations by uncertainties with the rock 
density/composition lies, at least, at a few percents level. 

All this have led us to a conclusion that we never know the 
muon energy losses better than with 1%-accuracy (actually, 
worse) and, consequently it makes no sense to waste efforts 
trying to reproduce energy losses with accuracy better than 
10~^ with an algorithm for the muon transport. So, there is 
a reasonable quantitative criterium for a muon propagation 
algorithm goodness: algorithm whose accuracy is better than 
1% can be accepted as a good one, otherwise it is bad. 



3 Choice of Vcut 



E„=100TeV: 



E„=10PeV: 




Fig. 5. Survival probabilities p vs Vcut obtained as a result of MC 

simulation for monoenergetic muon beams with initial energies Es 
= 1 TeV (1st column of plots), 10 TeV (2nd column), 100 TeV (3rd 
column) and 10 PeV (4th column). Each column contains six plots 
which correspond to six depths D. See text for more details. 



Among all simplifications that may be done when developing 
an algorithm for the muon propagation (and that affect the 
algorithm accuracy) one is obligatory. As number of muon 
interactions per unit of path is practically infinite it is impos- 
sible to simulate all acts and one has to set some threshold 
for relative energy transfer Vcut (v = i^E^/E^, where AE^ 
is energy which is passed by muon of energy E^^ to either 
real or virtual photon in a single interaction) above which 
muon energy losses are treated by the direct simulation of 
AE for each interaction ("stochastic" part) and below which 
energy losses are treated by means of stopping-power for- 
mula ("soft" part) that can be obtained by integration of dif- 
ferential cross-sections for the muon interactions: 



(A is a mass number of the target nucleus, Na is the Avo- 
gadro number). Number of interactions to be simulated per 
unit of muon path depends on Vcuu growing roughly as N^nt oc 
v~^^ along with computation time. So, it would be desirable 
to choose Vcut as large as possible. But, on the other hand, in- 
crease of Vcut affects the simulation accuracy since some part 
of statistical fluctuations in energy losses goes out of simu- 
lation and is not accounted for. Thus, the question is how 
large value ofvcut nray be chosen to keep result within 1%- 
accuracyl This problem was discussed in literature (Nau- 
movefa/., 1994), (Lipari and Stanev, 1991), (Antonioh 



al, 1997), (Lagutin et al, 1998) but, in our opinion, more 
careful analysis is lacking. 

Let's firstly consider propagation of monoenergetic muon 
beams. Fig. 5 shows results on propagation (survival proba- 
bilities p vs Vcut) for four beams with initial energies 1 TeV, 
10 TeV, 100 TeV and 10 PeV through slant depths 3.2 km, 12 
km, 23 km and 40 km, respectively, in pure water. Closed cir- 
cles represent results obtained with knock-on electrons (that 
appear as a result of ionization with large energy transferred 
to an atomic electron) included in direct simulation, open cir- 
cles stand for completely continuous ionization that is treated 
simply by Bethe-Bloch formula. Solid lines on each panel 
correspond to Vcut = lO^"' and all energy losses multiplied 
by factors 0.99 (upper line) and 1.01 (lower line), thus these 
lines limit "±l%-energy-losses-error-band". Horizontal dot- 
ted lines correspond to Vcut = 10^^ and cross section for ab- 
sorption of a real photon by a nucleon parametrized accord- 
ing to (Breitweg et al, 1999) instead of parameterization 
(Bezrukov and Bugaev, 1981a) that was mainly applied at 
presented simulations. 

One can conclude as follows: 

(a) In most cases with except for lower row and the left col- 
umn (that corresponds to low survival probabilities and 
low muon initial energies, respectively) survival proba- 
bilities he mainly within ±l%-band. 

(b) The difference between survival probabilities for two 
models of ionization - partially fluctuating (above Vcut) 
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and completely "soft" - is the less appreciable the larger 
muon energy is. 

(c) Parameterizations (Bezrukov and Bugaev, 1981a) and 
(Breitweg et ai, 1999) do not differ noticeable from 
each other in terms of survival probabilities. 

(d) For Vcut < 0.02-^0.05 there is almost no dependence 
of survival probabiUty on Vcut with except for very last 
part of muon path where survival probabihty becomes 
small. 

Now let's go to the more realistic case with atmospheric 
muons. Fig. 6 represents intensities of vertical atmospheric 
muon flux at 8 depths D of pure water as functions of Vcut 
as obtained by simulation with muons sampled according to 
sea level spectrum (Khmushin et al, 2001): 
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Meaning of closed and open circles, as well as solid and dot- 
ted lines is the same as for fig. 5. Dashed lines on panels for 
D < 5 km correspond to intensities which were computed 
for all energy losses treated as completely "continuous" (no 
fluctuations included). Dash-dotted lines show intensities of 
vertical muon flux simulated with muons sampled according 
to the Gaisser sea level spectrum (Gaisser, 1990): 
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Fig. 6. Intensities of vertical atmospheric muon flux at different 
depths of pure water vs Vcut as obtained by MC simulation. See 
text for more details. 
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General conclusions are qualitatively the same as for mono- 
energetic muon beams but quantitatively the influence of Vcut 
and model of ionization is much weaker. One can conclude 
the following: 

(a) Computed muon flux is strongly affected by accounting 
for fluctuations in energy losses: muon flux intensity 
computed by means of stopping-power formula is less 
comparing with MC simulated flux by «10 % at 3 km 
and by w 20% at 5 km. At the depth of 20 km muon 
flux computed with ignorance of fluctuations is only 10 
% of simulated flux. 

(b) 1% -uncertainty in muon cross sections plays the princi- 
pal role for resulting error (all simulated data he within 
±l%-band with no exceptions). This error has a ten- 
dency to grow with depth from ±2.5% at 1 km to ~ 

±15% at 20 km. 

(c) Difference between muon spectra (Klimushin et ai, 
2001) and (Gaisser, 1990) leads to uncertainty from 
-4% at 1 km to +16% at 20 km. 

(d) Error which appears due to simphfied, entirely "contin- 
uous" ionization lies, cormnonly, at the level of 2-^3%. 

(e) Dependence of simulated muon flux intensity upon Vcut 
is the most weak one comparing with other error sources. 



Thus, we can make the following conclusions: to repro- 
duce muon energy losses with an accuracy better than 10~^ 
it is quite enough to account only for fluctuations in energy 
losses with fraction of energy lost being as large as v > 
Vcut ~ 0.05-=-0.1 (commonly, values Vcut = 10~^-^10~^ are 
incorporated in majority of algorithms for the muon prop- 
agation by now (Lipari and Stanev, 1991), (Antonioli et 
al., 1997), (Rhode and Chirkin, 2001)). At least, in case 
with atmospheric muons only radiative losses may be simu- 
lated while ionization may be treated as an entirely continu- 
ous process without any lost of accuracy. In spite of remark- 
able dependency of simulated survival probabilities upon ac- 
counting or not accounting for fluctuations in ionization (see 
fig.5) it seems to be correctly to treat ionization by com- 
pletely continuous way also for neutrino-produced muons, as 

(a) differences in survival probabilities appear only for rela- 
tively low energies; (b) they are of few percents only; (c) they 
are of both signs and so, are expected to be partially self- 
compensated; (d) the uncertainties with i/ N cross-sections 
are still much larger, in any case (Agrawal et al., 1996), 
(Gaisser, 2001), (Lipari, 2001), (Naumov, 2001). With such 
parameters there are only several muon interactions to be 
simulated per 1 km w.e. 

We would like to emphasize that uncertainties with atmo- 
spheric muon flux intensity that are caused by only ±1% 
variations in energy losses are much higher than 1 % (see item 

(b) above). This is due to sharp power-law shape of surface 
atmospheric muon spectrum and due to the fact that source 
of atmospheric muons is far away from an underwater/ice 
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Fig. 7. Comparison between input models for the muon energy 
losses (lines) and simulated energy losses (markers) for muon trans- 
port codes PROPMU [version 2.01, 18/03/93], MUSIC [version 
for pure water with bremsstrahlung cross-sections by Kehlner- 
Kokoulin-Petrukhin, 04/1999], and MUM [version 1.4, 12/2001] 
in pure water. The energy range 25 GeV-^-l PeV is divided into two 
parts to keep linear scale: the upper plot stands for energy losses vs 
Eft = 25 GeV-^ 1 TeV, the lower one does for energy losses divided 
by energy vs = 1 TeV-^1 PeV. The model for PROPMU energy 
losses is taken from (Lohmann et al, 1985) as referenced in (Li- 
pari and Stanev, 1991). Curves for MUM's and MUSIC'S energy 
loss models are very close to each other and therefore are hardly 
distinguished. 

detector and their flux may only decrease when passing from 
the sea level down to detector depth. The source of muons 
that are produced in vN interactions is uniformly distributed 
over water/ice and/or rock both out- and inside the array. So, 
uncertainties with energy losses should affect the results for 
neutrino-generated muons weaker. Indeed, intensity of the 
muon flux which accompanies the neutrino flux in a medium 
is proportional to the muon range that, in turn, is inversely 
proportional to energy losses and, consequently, an error in 
energy losses leads to an equal error with an opposite sign in 
simulated flux of muons which are bom in u A'^-interactions. 



4 "Inner accuracy" of muon transportation algorithm 

Any muon MC propagation algorithm consists of a set of pro- 
cedures on numerical solution of equations, interpolation and 
integration. All these procedures are of finite accuracy and, 
consequently, the incoming model for muon interactions is 
somewhat corrupted by them. Having a set of formulae for 
the muon cross-sections for bremsstrahlung, e+e^-pair pro- 



Elg. 8. The same as in fig.7 with pure water, MUM [version 1.4, 
12/2001] wd 2 versions of PROPMU [versions 2.01, 18/03/93 and 
2.1, 01/2000] with different values of Vcut and different random 
generators. 

duction, photonuclear process, and knock-on electron pro- 
duction it is easy to obtain the formula for the total aver- 
aged muon energy losses by integration of differential cross- 
sections over all range of kinematically allowed v: 

Vmin 

But energy losses as they are simulated by an algorithm will 
differ from (4) due to "corrupting" application of numeri- 
cal procedures. Two algorithms with the same formulae for 
cross-sections at the input will reproduce the muon energy 
losses differently. The difference between simulated and cal- 
culated energy losses collects all errors that are contributed 
by each step of algorithm and thus, is a good quantitative 
criterium for its inner accuracy, whose contribution to the 
resulting error must not exceed 1%, as it was shown in the 
Sec. 2. 

The trivial test may be suggested to check inner accuracy 
of any muon transport algorithm: let's simulate propagation 
of n muons (we used n =10^ for test reported below) with 
energy Eq over a short distance Al "measuring" the final en- 
ergy El for each i—th muon at the end. Then, simulated 
energy loss may be computed as 

1=1 

Comparing (5) with (4) it is possible to determine a good- 
ness of the given algorithm applying, for example, criterium 
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Fig. 9. The same as in fig.7 and fig.8 with standard rock (p = 2.65 
g/cm^, A = 22, and Z = 11), MUM [version 1.4, 12/2001] and 
PROPMU [version 2.1, 01/2000]. 



abs{{5) - (4)) < 1%. 

Fig.7, fig.8, and fig.9 show results of such a comparison 
for three muon transport codes: PROPMU (Lipari and Stanev, 
1991), MUSIC (Antoniohef a/., 1997), andMUM (Sokal- 
ski et al. , 2001). Results for pure water are presented in fig.7. 
In case with MUM the difference abs{{5) — (4)) never ex- 
ceeds 1%, for MUSIC it reaches 2-^3%. Inaccuracy of the 
PROPMU code is much worse and in a range < 10 TeV 
is up to 25%. We try 2 versions of PROPMU with differ- 
ent values of Vcut and different random generators (fig.8). In 
some cases iimer accuracy becomes better but, in any case it 
stays worse than 10% that is unacceptable according to con- 
clusion of Sec. 2. Fig.9 reports the results of accuracy test for 
MUM and PROPMU with standard rock. QuaUtatively, these 
results are the same. 

In fig. 10 differential spectra for atmospheric muons at dif- 
ferent depths in pure water are presented as simulated with 
PROPMU, MUSIC and MUM. Muons at the sea level were 
sampled according to the same spectrum (Klimushin et al., 
2001). Also parameterizations (Klimushin etal., 2001) and 
(Okada, 1994) are shown. Simulations with MUSIC and 
MUM give practically the same results because both codes 
reproduce almost the same energy losses (see fig.7). MU- 
SIC'S and MUM's spectra coincide with parameterization 
(Klimushin et al, 2001) that is based on the same sea level 
muon spectrum as was used for simulation presented on the 
plot and on muon propagation with MUM. Parameterization 
(Okada, 1994) is lower than KBS, MUM and MUSIC (up to 
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Fig. 10. Vertical muon depth spectra in pure water at depths 
D = 1 km, 3 km, 6 km and 10 km, as simulated with PROPMU, 
MUSIC and MUM for atmospheric muons sampled according to 
sea-level spectrum (Klimushin et al., 2001). For a comparison, pa- 
rameterizations (Klimushin e^a/., 2001) (KBS) and (Okada, 1994) 
(Okada) are given, as well. On panel for D = 1 km the spectrum 
which resulted from simulation with MUM with all cross-sections 
multipUed by 0.9 is also presented (MUM*, open circles). 



18% in terms of integral muon flux at D = 1 km) at relatively 
shallow depths and becomes higher starting with I? ~ 5 km. 
The reason is that it is based on very hard sea-level muon 
spectrum with index 7 = 2.57 (Miyake, 1973), while KBS 
parameterization adopts a spectrum with 7 = 2.72. As a re- 
sult, there are less muons with relatively low energies (which 
contribute the main bulk to the muon flux at low depths) and 
more muons of high energies which contribute to the muon 
flux at larger depths. The reduced muon energy losses in 
PROPMU result in the depth muon spectra which (a) are sig- 
nificantly higher (31%, 30%, 27% and 17% in terms of inte- 
gral muon flux at the depths D = I km, 3 km, 6 km and 
10 km, correspondingly) and (b) are expanded to the low 
energies. The spectrum for D = 1 km was simulated also 
with "corrupted" version of MUM: all muon cross-sections 
along with Bethe-Bloch ionization formula were multiplied 
by a factor 0.9, thus energy losses were reduced by 10% and 
became close to ones that are reproduced by PROPMU ^. 
The result is given on upper left panel of fig. 10 by open cir- 
cles (MUM*). A good agreement between PROPMU and 
MUM* represents a good quantitative cross-check for results 

'The majority of muons which compose the spectrum at 1 km 
depth have the energy in a range E < I TeV at the sea level. For 
this range PROPMU reduces the simulated energy losses by ~ 10% 
as can be seen in upper panel of fig.7. 
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on PROPMU inner accuracy. 



5 Conclusions 

The optimum value for inner accuracy of a MC code for the 
muon transport is close to 1%. In any case we never know 
the muon energy losses better. 

To keep such irmer accuracy it is quite enough to simulate 
the muon interactions with a fractions of energy lost that are 
as large as v > Vcut = 0.05-^0.1, treating the rest of energy 
losses (soft part) by means of stopping-power formula. 

The simple test may be applied to measure the accuracy 
of an algorithm for the muon transport by comparison of re- 
produced energy losses and incoming model for muon cross- 
sections. 

Such a test has been applied to three muon transportation 
codes. According to "1% goodness criterium" codes MUM 
and MUSIC can be accepted as accurate enough. Accuracy 
of the PROPMU algorithm is insufficient (down to 25%) and 
depends on medium, version of code and simulation parame- 
ters. It leads, for example, to Ri30% overestimation of atmo- 
spheric muon fluxes at depths l-j-5 km w.e. 
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